Magnetic field-induced soft mode in spin-gapped high-T c superconductors 
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We present an explanation of the dynamical in-gap spin mode in L^-^Sr^CuCU (LSCO) induced 
by an applied magnetic field H as recently observed by J. Chang et aL— Our model consists of a 
phenomenological spin-only Hamiltonian, and the softening of the spin mode is caused by vortex 
pinning of dynamical stripe fluctuations which we model by a local ordering of the exchange inter- 
actions. The spin gap vanishes experimentally around H = 7T which in our scenario corresponds 
to the field required for overlapping vortex regions. 
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The cuprate superconductors (SC) arise from doping 
an antiferromagnetic (AF) Mott insulator. At half-filling 
the spin spectrum of La2Cu04 is quantitatively described 
by a spin- 1/2 Heisenberg models At finite doping, how- 
ever, the nature of the magnetic fluctuations and their 
importance for SC remains controversial. At present, 
the so-called "hour-glass" dispersion observed in inelas- 
tic neutron response appears universal whereas details 
of the low-energy spin fluctuations vary between the 
compounds^ In the optimal and overdoped regimes an 
interplay between magnetism and SC has been revealed 
by the opening of a spin gap which scales with the SC 
transition temperature T c . This is in contrast to the 
underdoped regime where LSCO and Bi2Sr2CaCu2 0s+5 
(BSCCO) are know to exhibit spin freezing well into the 
SC domefi&&! whereas YBa 2 Cu 3 6+x (YBCO) reveals 
a static signal only at very low doping^ 

In the quest of obtaining a better understanding of 
these materials, the effect of an applied magnetic field 
H has been extensively used, also in neutron scatter- 
ing experiments. In the underdoped regime a magnetic 
field applied perpendicular to the Cu02 planes enhances 
static incommensurate (IC) stripe order, which exist at 
H = presumably due to impurities ! 10 ' 11 This enhanced 
signal at the IC positions q* [quartet of peaks near (ir, tt)] 
has been seen both in underdoped L^-zSr^CuO^ ' 12 ^ 3 
La2Cu04 +y r^ and very recently also in YBa2Cu30g.45i£ 
A similar enhancement of stripe order can be obtained 
at H — by the use of impurity substitution. 15 : 16 ' 17 : 18 
Furthermore, experiments have shown that spin-gapped 
samples can transition from a SC phase to a coexisting 
SC and IC stripe ordered phase by use of a magnetic 

In LSCO, static order is absent for doping levels be- 
yond approximately x ~ 0.13/ 4 - The inelastic spin exci- 
tations are, however, still characterized at low energy by 
the same IC wavevectors but a spin gap of order ~ 3 — 7 
meV develops at T < T C 2&2L2Z22. The iJ-dependence of 
the low-energy inelastic neutron response has also been 
subject of intense experimental investigations. Lake et 
al£^ reported a softening of the spin mode in LSCO 
(x = 0.163) revealed by an in-gap mode observed with 



H = 7.5T compared to a fully spin-gapped spectrum 
for H = 0. Similar results have been obtained at larger 
doping levels i 25 : 26 More recently detailed inelastic neu- 
tron scattering experiments studied the iJ-dependence 
of the magnetic spin gap in slightly underdoped LSCO 
(x = 0.145)i At this doping, a critical field of H c = 7T 
is required to tune the system from a SC state into a 
phase with coexisting SC and long-range IC stripe order. 
At applied fields < H < 7T the spin gap is diminished 
and an in-gap spin mode is observedJ^l This transition 
is very reminiscent of the dynamic neutron response seen 
e.g. by Kimura et aL 16 in Zn-doped LSCO (x = 0.15). 

The presence of SC regions appear important for the 
existence of enhanced magnetic response at these rela- 
tively low applied fields H ~ — 10T; only when vortices 
can act as additional pinning centers for nucleation of 
stripe regions do they lead to an enhanced magnetic re- 
sponse. This agrees qualitatively with: 1) the fact that 
the enhanced signal is seen at T < T c , and 2) an ab- 
sent [a negligible] magnetic field effect in non-SC [weakly 
SC] samples! 27 ' 28 : 29 It is not necessary for the vortices 
to form an ordered lattice which also appear absent at 
x = 1/8 in LSC03 Theoretically, the existence of AF 
order induced by vortices was first discussed within the 
SO (5) theory of the cuprates i 30 : 31 ' 32 Later several mod- 
els studied how vortices may nucleate magnetic regions 
due to a general competition between SC and stripe 
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One way to model the stripe phase is in terms of cou- 
pled spin ladders! 40 ' 41 : 42 : 43 : 44 ' 45 In such spin-only models, 
the charge carriers are assumed important only for renor- 
malizing the exchange couplings between localized spinful 
regions, and the Hamiltonian is of the Heisenberg form. 
Clearly this approach is phenomenological, and should be 
considered an approximate effective model. Nevertheless, 
this approach has been very successful is describing e.g. 
the universal "hour-glass" magnetic dispersion. 

Here, motivated by the recent experimental findings of 
Chang et air- we study theoretically the effect of a mag- 
netic field on the low-energy (gapped) spin fluctuations. 
By extending the coupled spin ladder approach to include 
the effect of vortices, we find a field-induced mode inside 
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the gap similar to experiment. The vortices are modeled 
by local regions of exchange couplings J v different from 
the bulk J;,. Thus, the effective Hamiltonian reads 

H = JbSi ■ Sj + J„Si • Sj, (1) 

b(ij) v(ij) 

where (ij) denote neighboring spins, and b and v re- 
fer to effective exchange couplings far away from (in 
the vicinity of) the vortices J b = { J^ a , j£ lter } ( J v = 
{j; ntra , J™ ter }). For simplicity, we assume that the intra- 
ladder exchange couplings are unaffected by the vortices 
jmtra _ jintra _ anc j fa&t only the inter-ladder cou- 
plings may differ between b bulk and v vortex regions. 

Self-consistent microscopic studies of the Hubbard 
model have shown that, depending on input parameters, 
impurities and vortices can induce uni-directional stripes, 
approximately rotational invariant stripes (ARI stripes), 
or checkerboard patterns in the surrounding magnetiza- 
tion and hole density 35 i 36 ' 46 i 47 Furthermore, the induced 
spin density is found to be modulated with a period close 
to eight lattice spacings i 35 i 36 Consistent with these re- 
sults we model the vortex regions by ARI stripes and 
arrive at a picture of the spins near the vortices as shown 
in Fig. [TJ Clearly, this spin arrangement is idealized, 
but we believe it gives a reasonable description of the 
structure near the vortices. For example, as shown in 
Ref. [H, stripe arrangements similar to Fig. [T] agree with 
the checkerboard local density of states modulations ob- 
served near vortices in BSCC04& We stress, however, 
that the particular choice of induced order shown in 
Fig. [1] is not important for the following discussion; sim- 
ilar results can be obtained from uni-directional stripes. 
Checkerboards, on the other hand, do not appear to 
be consistent with experiment^ The important point is 
that vortices in underdoped cuprates nucleate stripe or- 
der in a surrounding halo characterized by a new length 
scale £jc> which lies between the SC coherence length £ 
and penetration depth A for these materials, £ < £jc < A. 
Typical spatial extent of the vortex regions used in our 
simulations is taken to be ~ 50 x 50 lattice sites consistent 
with estimates of the correlation length £xc ~ lOOAji^ 

The value of the inter-ladder exchange coupling J™ ter 
is restricted by the fact that for high magnetic fields 
the vortex regions overlap and long-range IC order is 
observed^ implying that | J™ tcr | should exceed the criti- 
cal value J™ ter « —0.4 J (for 2-leg ladders^ 4 -) for appear- 
ance of long-range order. Thus, in the following we sim- 
ply fix 4 ntcr = -0.6 J, i.e. well below J™ tel . On the other 
hand, the absence of long-range spin order at H = im- 
plies that the bulk inter-ladder exchange couplings J™ ter 
are weaker than this critical value. While it is possible 
to obtain the results reported here by assigning a small 
value to j J™ ter | , we prefer instead to view J™ ter as an ef- 
fective coupling obtained by randomly diluting bonds of 
strength J^ ntcr . Thus, some of the J™ tcr are while oth- 
ers are —0.6 J, mimicking the disordered nature of stripes 
in the bulk & i 10 i 49 i 50 i 51 ' 52 i 53 By contrast, we assume that 
the primary effect of a magnetic field is to introduce vor- 
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FIG. 1: Arrangement of stripes around a vortex core. The ar- 
rows indicate spins and show an ordered configuration accord- 
ing to the signs of the couplings. The dots are nonmagnetic 
sites. The AF interaction between neighboring spins on the 
same ring is J^ ntra while spins on different rings are coupled 
by a ferromagnetic interaction J™ tel which are indicated, for 
a few spins, as straight lines in the lower left corner. 

tex regions of ordered stripes where all the inter-ladder 
bonds have strength J™ tor = — 0.6J. The motivation 
for this effective " exchange-ordering" nature of H comes 
from recent quantum oscillation measurements which in- 
dicate that an external magnetic field cause exactly this 
kind of ordering of the stripes resulting in a severe Fermi 
surface reconstruction.— 

In order to find the spin excitation spectra we simu- 
lated Heisenberg spin systems as depicted in Fig. Q] at a 
temperature T = 0.01J using the quantum Monte Carlo 
stochastic series expansion technique with directed-loop 
updates. 55 This technique yields high quality results for 
the imaginary-time correlation function which is con- 
tinued to real frequencies using the Average Spectrum 
Methodp& where the average over all possible spectra is 
taken weighted by how well each spectrum fits the data. 
This approach performs at least as well as MaxEnt meth- 
ods for high quality imaginary-time data^i 

The bulk part Sb(q*,uj) and the vortex part S v (q*,oj) 
of the structure factor were simulated separately. For 
the bulk we diluted the arrangement of J^ ntcr bonds ran- 
domly prior to performing the simulations. The structure 
factor Sb(q*,u>) was extracted for each disorder realiza- 
tion and the average was taken. With an interladder 
bond dilution fraction of 30% we found that Sb(q* , u>) has 
a sharp peak around uj = 0.1J as well as a broader high 
energy peak at u> = 0.2 J as seen from the inset in Fig. [2] 
This two-peak structure at the IC position q* is con- 
sistent with inelastic neutron scattering measurements 
on spin-gapped LSCO* 5 ^ Without a microscopic model 
including realistic disorder concentrations it is hard to 
estimate the correct degree of bond disorder within the 
present spin-only approach. However, the position of the 
low energy peak depends on the amount of bond dilution, 
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FIG. 2: Magnetic structure factor Si(q*, u>) at the IC wavevec- 
tor q* = (37r/4, n) as a function of energy lo. This spin 
spectrum is obtained from a disordered set of coupled 2-leg 
spin ladders mimicking the neutron response from disordered 
stripes at H = 0. The inset shows the same quantity for a 
larger energy range. 
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FIG. 3: Total magnetic structure factor S(q*,Lo) at the IC 
wavevector q* = (37r/4, 7r) as a function of energy lo including 
the effects of an applied magnetic field. The inset shows the 
same quantity for a larger energy range. 

less dilution moves the peak down in energy, and we have 
simply chosen 30% in order to reproduce a spin gap of 
order 10% of J as reported in Ref. QJ. 

For the vortex regions we use the non-diluted configu- 
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ration to obtain S v (q*,Lij) which has a structure similar 
to Fig. [5] except that the low-energy peak has moved fur- 
ther down in energy. For J™ tcr = —0.6 J the lowest peak 
is roughly at lo — 0.02J (see Fig. [3]). This corresponds 
to lo ~ 2meV in agreement with Ref. 1 using recent es- 
timates for the exchange coupling J ~ lOOmeV in LSCO 
(x = 0.14)i££ The peak position is dependent on both 
the finite size of the vortex region and the value of J™ ter ; 
as expected smaller vortex regions (or smaller value of 
|J™ ter |) moves the peak upwards in energy. Taking into 
account that for H = 2.5T (as used in Ref. [l]) the bulk 
region contributes roughly 15 times more to the total 
S(q*,uj) than the vortex regions, we compose the total 
structure factor S{q*,Lo) = (15S b {q*,Lo) + S v {q*, w))/16 
shown in Fig. [3J As seen, H induces an in-gap mode 
similar to the experimental results of Refs. 0] and [24]. 

At low magnetic fields [H < 3T in LSCO (x = 0.145)], 
an increase of H enhances the vortex density, but the size 
of each vortex region is presumably unchanged resulting 
in a fixed energy of the in-gap spin mode. In this field 
range the spectral weight of the in-gap mode should in- 
crease as a consequence of an increased weighting of the 
total structure factor by the denser amount of vortex re- 
gions. By contrast when vortex regions overlap, the mode 
will rapidly move to zero energy becoming a true Bragg 
signal. Vortex regions of size ~ 100A will start overlap- 
ping at H ~ 6T in agreement with the closing of the spin 
gap found in Ref. [0 near this field strength. Theoreti- 
cally, the generation of a Bragg peak happens if J™ ter is 
larger than the critical value J™ ter needed for long-range 
stripe order. By contrast, if | J^. nter | < |</c ntcr |> a finite 
energy peak can still show up at finite field, but it will 
not move to zero as the vortex regions start overlapping. 

In conclusion we have proposed an explanation for the 
soft magnetic mode observed inside the spin gap of LSCO 
in terms of magnetic field-induced stripe ordered vortex 
regions. At zero applied magnetic field H = the Cu02 
planes are modeled as disordered coupled spin ladders 
known to reproduce the magnetic " hour-glass" dispersion 
seen by experiments. At finite field H ^ vortices are 
simulated by regions of ordered exchange couplings. Our 
calculation shows that stripe pinning by vortices gives a 
consistent picture of the in-gap spin mode observed in 
recent inelastic neutron scattering measurements. 
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